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We note that a simple two parameter description of lepton mixing is possible which 
reproduces the features that apparently emerge from global fits at the lcr level: if 
U e 3 is non-zero it implies that the solar neutrino mixing parameter sin 2 #12 is less 
than I by order |i7 e 3| 2 - If the CP phase 5 is around tt it implies that the atmospheric 
neutrino mixing parameter sin 2 #23 is less than ^ by order | C/ e 3 1 . The mixing scheme 
can be described by a 23-rotation appearing to the right of a tri-bimaximal mixing 
matrix. We quantify the excellent agreement of the scheme with data statistically, 
and comment on model building aspects. 
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INTRODUCTION 



All three mixing angles of the lepton sector are now known. The last step towards 
this marvelous achievement came from reactor neutrino experiments Double Chooz [l], 
Daya Bay Q| and RENO j^J. Combining the reactor data with other experiments ruled 
out vanishing U e3 at more than 7a C.L. [4-7]. At the Neutrino 2012 conference in 



June 2012, Double Chooz have presented new data with 3.1<r evidence for non-zero 
U e3 , and also Daya Bay has increased its significance to more than 7a, see the URL 



http : //kds . kek . jp/conf erenceTimeTable . py?conf Id=9151#20120604 . detailed for the 



slides. Moreover, additional data from T2K and MINOS can be included, and one can fit 
the overall data to the parameters in the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) lepton 
mixing matrix (ignoring possible Majorana phases) 



f C12C13 S12C13 s 13 e~ l5 \ 

U = -S12C23 - ci 2 S23Si 3 e 1<5 C12C23 - si 2 S23Si 3 e 1<5 S23C13 

\ S12S23 - Ci 2 C23Sl 3 e 1 ' 5 -C12S23 - Si2C2 3 Si 3 e 1 ' 5 C23C13 J 



In this short note we wish to give a possible interpretation of the emerging features of global 
fits. We will focus on the results from Fogli et al. f| , which spotlight the following interesting 
properties: 

a) \U e3 \ ~ 0.16 is sizable; 

b) solar neutrino mixing is described by sin 2 12 ~ | — 0(\U e3 \ 2 ), i.e. slightly less than |; 

c) atmospheric neutrino mixing is described by sin 2 6 23 ~ \ - 0(\U e3 \), i.e. significantly 
less than |; 

d) the CP phase 5 is around n. 

The precise fit parameters are quoted in Table I. Not all features are present in the results of 
the other groups 0,0,0] (which have not yet updated their results, and have partly different 
treatment in their atmospheric codes or do not even fit atmospheric data), and at the 2a 
level points b) and d) are absent. While there is no doubt about the value of \U e3 \, less-than- 
maximal atmospheric mixing and less-than-| solar mixing seem to be common features, at 
least at the la level. 

If these properties of lepton mixing survive the test of time, an interpretation in terms of 
a mixing scheme will without doubt be useful. We note in this work that there is a mixing 
scheme that can reproduce features a) - d). It has only two free parameters that can be 
adjusted to the observables, and possesses the following properties: 

i) if \U e3 \ is non-zero, sin 2 6*12 — | — 0(|£/ e 3| 2 ) is implied. This links features a) and b); 

ii) if 5 is around ir, sin 2 6 23 — \ — 0{\U e3 \) is implied. This links features c) and d). 
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TABLE I: Best-fit and estimated la values of the neutrino mixing parameters from Ref. |f| for the 
normal mass ordering (7723 > m 2 > mi, NH) and the inverted mass ordering (771,3 < m i < m 2, IH). 



Parameter 


sin 2 9i2 


sin 2 023 


sin 2 13 


5/tt 


NH 
IH 


o w +0 - 018 

U - JU ' -0.016 


0.38618EJ 

q oo 2 +0.039 

u.oyz_ 022 


0.0241 ±0.0025 
0244+ ' 0023 

U.UZ^_ Q 0025 


1.08%* 
1.091°;! 



The defining property for the PMNS matrix U is 



1 <r\ # #^ 



v 



(2) 



, # # 

where the elements in the second and third column can be obtained by unitarity. The 
entries of this matrix in the first column are exactly as in tri-bimaximal mixing (TBM) [s| , 
and hence this mixing matrix can also be obtained by multiplying the tri-bimaximal mixing 
matrix with a 23- rotation from the right. The phenomenology of the mixing scheme has 
been analyzed first in {^J (see also 1Q|), and is a variant of the so-called trimaximal mixing 
scheme, which is defined through a mixing matrix with only the second row as for TBM 
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.2/3 ^ 1(1 



\u, 



e3 I 



The 



13]. This scheme is however disfavored due to its prediction sin 612 - 
scheme that we propose to describe the current data has been written down first in 13], 
and in the convention of ^\ is called TMi. We revisit this mixing scheme in this short note, 
emphasizing that it is able to perfectly accommodate the above features a) - d), which seem 
to emerge from global fits. In Section [TT] we study its phenomenology, and in Section III II we 
discuss aspects of a possible theoretical background. 



II. PHENOMENOLOGY 



Consider the tri-bimaximal mixing matrix multiplied with a 23- rotation from the right: 

/ ./I ./I \ 



U = U TBM R2 3 (6^) 



The observables are in this case 

1 



5 / 



\u e3 \ 2 = 

sin 2 0i2 



sin 2 9 , sin 2 9 



23 



3 - sin 2 9 



1 

~ 2 
Jcp 



f 1 \ 

cos 9 e~ 1 ^ sin 9 
\0 -e^sinfl cos 9 j 



I sin 29 cosi/j 

3 - sin 2 9 ' 
sin 29 sin ip , 



6V& 



(3) 



(4) 



where the Jarlskog invariant is defined as Jcp = Im (UeiU^U^U*^ = Si 2 Ci2S 2 3C23C 13 Si3 sin 5. 
The main feature of the above mixing matrix is that the first column of U keeps the same 
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form as for TBM: 



/ Well 2 \ 


/2/3\ 


Ki\ 2 - 


= 1/6 


VIM 2 J 





(5) 



From these relations, as well as from Eq. (j3J), the observable mixing parameters and their 
correlation can be obtained. In particular, a consequence of Eq. 
Indeed, from \U e i\ 2 



is that sin 2 6 12 < 



3" 



| one finds 



sin 2 # 



12 



1-3 


\u e3 


12 


1 - | 


u e3 \ 


2 



(l-2|f/ e3 | 2 ). 



(6) 



The desired feature b), sin 2 9 12 = \ — 0(\U e3 \ 2 ), is reproduced. In Fig. [TJ we show the 



correlation between \U e3 \ and sin 2 #12, using the la range of the global fit results [6|. For 
simplicity, we use the results for the normal ordering, the difference to the results for the 
inverted ordering is insignificant. One finds that sin 2 # 12 lies between 0.315 and 0.318. 



The second independent condition in Eq. ([5]) involving \U t 



cos 5 tan 26 



1-5\U, 



e3 1 



23 



2V2\U e3 \ 
This relation can be written as 



l-3\U e3 \ 2 2V2\U t 



c3 I 



1/6 gives 

i-I 

2 



\U, 



1 3 



sin 2 



23 



1 1 lrr , 

- + — F= ^3 
2 



4 1 eci| 



cos 5 . 



(7) 



(8) 



Hence, if cos 5 < and |cos5| = 0(1), the desired features c) and d) are reproduced: 
sin 2 6*23 = \ — 0{\U e3 \). We plot in Fig. [TJ the correlation between sin 2 6 23 and 5. To 
reproduce the la range of 8 23 , 6 should lie between 1.307T and 1.387T. It may be useful to 
express the Dirac CP phase in terms of the parameters 9 and ip as 

^/2{b + cos26) sin?/; 



sin 6 



(9) 



'39 + 20 cos 26 + 13 cos A6 - 24 cos 2ifj sin 2 2d 

We also display in Fig. [TJ the correlation between sin 2 #23 and sin 2 # 12 , when \U e3 \ is varied in 
its la range. 



For further illustration on how nicely the mixing scheme describes current data, we per- 
form a somewhat naive statistical fit. To this end, we take the two parameters 9 and ip as 
independent, and compare the TMi predictions to the experimental data with a x 2 function 

x ^j : iPi_ff ] (10) 

i 1 

where p\ represents the data of the z-th experimental observable, <Tj the corresponding la 
absolute error, and pi the prediction of the model. The experimental values of the neutrino 
mixing angles are taken from Table I, which as mentioned above include the data presented at 
Neutrino 2012. Note that the global-fit data slightly differ for normal and inverted neutrino 
mass orderings. 
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FIG. 1: Correlation of observables for our mixing scheme: the upper plot shows sin 2 Q\2 vs. |C/ e 3|, 
the middle plot sin 2 #23 vs. 5 and the lower plot sin 2 623 vs. sin 2 612 for different values of 5. Except 
for the lower plot, where the displayed range of 623 is larger, the plots cover the allowed la ranges 
of the parameters in case of a normal mass ordering [f|. The black solid lines are the best-fit points, 
the black dashed line in the lower plot the la range. 
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FIG. 2: The allowed region of the TMi parameters 9 and tp at la, 2a, and 3a C.L. for the normal 
mass ordering (left) and inverted mass ordering (right). The black asterisks denote the best-fit 
values. 



In Fig. El we present the allowed regions of 9 and ip at la, 2a, and 3a C.L., denned as 
contours in A% 2 for two degrees of freedom with respect to the x 2 minimum (Xmm — 1-1 
for NH and xLn - °- 9 for IH )- Note that these 

X 2 minima show that the scheme describes 
the data excellently. Since \U e %\ is firmly connected to 9 in TMi, it is restricted to a narrow 
range between 0.077T and O.Itt, and the best-fit value of 9 lies close to n/12. The constraints 
on i/j are not as strong as for 9 due to the less precise determination of 5 and #23 . The best-fit 
values of 9 and ip are (0.087vr, 0.32tt) for the NH case, and (0.087tt, 0.32tt) for the NH case. 

We further show in Fig. [3] the predictions for the mixing angles and the Dirac CP phase 
S for the NH case (upper panel) and the IH case (lower panel). Again, the data are very 
well reproduced. The best-fit values of sin 2 #12 and sin 2 6*23 are found to be (0.317, 0.382) 
for the NH case, and (0.317, 0.385) for the IH case. The right column shows that the TMi 
prediction on 6* 13 is in good agreement with experiments, and that 5 tends to lie at the upper 
end of its allowed range. We find the best-fit values of sin 2 #13 and 5 to be (0.0242, 1.34-71") 
for NH and (0.0245, 1.3571") for IH. For comparison, we also show in dotted contours the 
la parameter ranges without considering the experimental data on 5. In such a case, the 
allowed parameter spaces are symmetric with respect to 5 = n, and is confined to be either 
around = 1.37T or 5 = 0.7ir, implying the predictive power of the scenario. All in all, three 
mixing angles together with one CP phase are all compatible with experimental data within 
la C.L., though there are only two free parameters in TMi. 
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FIG. 3: The allowed region of the physical observables at la, 2a, and 3a C.L. for the normal 
mass ordering (upper panel) and inverted mass ordering (lower panel). The shaded areas are not 
permitted by the model. For comparison, the allowed la ranges |f| of the parameters are also 
indicated by using green vertical and horizontal bars. The dotted contours are the la ranges 
without considering the experimental data on S. 

III. THEORY 



One may ask a question arising from the previous analysis: what is the underlying symmetry 
behind TMx? In this section we wish to give some comments on the possible theoretical 
background and in particular we will show two examples in the framework of and £4 
flavor symmetries. 



A. Symmetry behind the TMi 



The original trimaximal mixing, defined by (\U P ?.\ 2 , lET^I 2 , \U^\ 2 ) T = (1/3, 1/3, 1/3) T , has 
been discussed in the framework of flavor symmetry models 10, ll2H2l|. Models that can 



s 



lead to the mixing scheme dealt with in this paper are, to the best of our knowledge, discussed 
only in 13|, |22| . The mass matrix that reproduces our mixing scheme is: 



/ A 



m,. 



U* W 



B + C 



B-C 



\ 



■ \{A + B + D + AC) \{A + B-D) 
\ ■ ■ \(A + B + D-AC) ) 



foe -c \ 

■ 2C 

where we identify (with c = cos 9 and s = sin 9) 



11] 



A = | (2 mi + m 2 Cg e 
B = \ {-m x + m 2 c 2 e~ 2ia + m 3 s 2 e e 2i ^~^) , D = m 2 e ~ 2i ^ + ^ s 2 + m 3 cj e~ 2if} 



+ m 3 s 2 e e 2i ^~^) , C = (m 2 e- i W' +2Q ) - m 3 e^- 2 ^) s e c e 



Here we have used the form of the PMNS matrix in Eq. ([3]), including for completeness the 
two Majorana phases. We see that the yU-r symmetry is broken by the extra terms involving 
C. In the limit of C = we would have tri-bimaximal mixing. Namely, 



m 



Hc=o 



m 



TBM 



C/tbm diag (A - B, A + 25, D) U, 



T 

TBM 



(12) 



Note that C is proportional to sin 9, and thus expected to be somewhat smaller than A, B, D. 

The eigenvalue A — B has an eigenvector ^g(2, —1, — 1) T , corresponding to the invariant 
column in TMi. This symmetry of the neutrino mass matrix can be characterized by a 
unitary matrix Gi, satisfying the relation m u — > Gjm v Gi, in which G\ defines a Z 2 group, 
i.e. G\ = 1, and it is given by 



/ 1 -2 
-2 -2 1 
-2 1 -2 



2\ 



(13) 



The third column of the mixing matrix allows us to define another Z 2 symmetry, under 
which the mass matrix is invariant, 



G 2 



/ 2 + c 2B 

s g (y/6t 



e**ce - 2s 6 



2s e ) 
-2s e ) 



Sg(2y/6c e COS^ + S e ) 



e^c e 



2s e ) \ 



\(l + 5c 26 i - 2zA/6sin?/;s 2 6i) 
[1 + 5c 2 e + 2i\/Q simps 2 e) — s e (2\/Qc e cos ip — s$) J 



Therefore, the neutrino mass matrix is invariant under the Z 2 x Z 2 transformation generated 
by Gi x G 2 - Note that we are working in the basis in which the charged lepton mass matrix 
is diagonal and non-degenerate. Therefore, any unitary matrix F that commutes with the 
charged lepton mass matrix must be diagonal with unit moduli in all its entries. For the three 
generation case, there are three distinct F: diag(l, to, u 2 ), diag(w, 1, uj 2 ) and diag(w 2 , 1, u) 
with u = exp(27ri/3). This allows one to construct a flavor group Q generated by F, G\ 



and G 2 , i.e. Q = {F, G\, G 2 }. The generator G\ is well motivated in severa 
groups, including A4, S4, and other higher order discrete groups (see Ref. 



flavor symmetry 



13| for a detailed 
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discussion, where in fact S4 was proposed for the scheme under consideration). In contrast, 
Gi cannot be simply embedded into small discrete groups due to the rotation angle 9. For 
6 = the matrix G 2 becomes the generator of ji-r symmetry, 



G 2 {6 = 0) 



/ 1 \ 
1 
1 



(14) 



For 9 = 7r/12 (which is close to the best-fit point) and <p = ir, one has 



G 2 {e = tt/12) 



/ I(4 + V3) ^(-4 + 2^-^) n (-4 + 2^/3 + VQ)\ 
(-4 + 2^-^) ^(2-^3-2^6) 
4 + 2^+^) ^(2 + 5^) 



12 



12 



(2 + 5^3) 



1 

4^ 



,(15) 



V6 



which actually gives \U e3 \ 2 ~ 0.025, sin 2 12 ~ 0.318 and sin 2 ~ 0.709. Both 13 and 9 12 
are well within the current la ranges, while ^23 deviates from its best-fit value at more than 
2a C.L., which might be improved once an explicit model with perturbations is constructed. 



B. Realization in flavor symmetry models 

As we mentioned above, the TMi mixing scheme can be viewed as a modification to the 
TBM mixing pattern by multiplying a 23-rotation from the right to £/tbm- In this sense, 
the rotation matrix R2z{0,ip) could also be viewed as a perturbation to the exact TBM 
mixing pattern. Alternatively, and this is what happens in the two short examples we are 
about to give, one can note that the mass matrix for TM4 is the TBM mass matrix plus 
a simple additional term, see Eq. (Hip . It is therefore possible that we modify a successful 
model leading to TBM, and add additional flavons and particles to it which give precisely 
the required form for TMi. 



In the original Altarelli-Feruglio model 23|, |24f , the left-handed lepton doublets are as 



signed to a three dimensional representation 3 under the tetrahedral group A4, whereas the 
right-handed lepton fields transform as 1, 1" and 1', respectively. In addition, three sets 
of flavon fields <px, <ps an d £, transforming as 3, 3 and 1 under A4, are also introduced 
together with the vacuum expectation values 1 (VEVs): (<pt) = (vt, 0, 0), (<ps) = ( v s, vs, vs) 
and (£) = u, respectively. At leading order, the A4 invariant Lagrangian contains terms like 

£ = Ve^e c (ip T £) h d + y^fi c (ip T £)' h d + Vr^r c (ip T i)" h d 
+7^^aC {Zhjhu) + -£jx b (p s (ih u £h u ) + h.c, 

where y a and Xi denote the corresponding Yukawa couplings, A is the cut-off scale of the 
theory, and two Higgs doublets h u and h d with VEVs v u and v d are assumed to be invariant 



1 Note that we assume the solution of the vacuum alignment could be achieved. In fact, several methods 
(e.g. by using driving fields) have been proposed to explain the different alignments, and these will not be 
discussed here. 
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under A 4 . Note that there is also an additional Z s symmetry in the model, which decouples 
the charged lepton and neutrino sectors 



23 



241 ] . By inserting the VEVs of the flavon fields, 



one finds that the charged lepton mass matrix is diagonal at leading order, i.e. 



m 



(o) 



VdVT 



1 A 

whereas the neutrino mass matrix is given by 



diag (y e , y T ) , 



(17) 



with a = 
Utbm as 



2x Q ^ and b 



m 



ZX *> A ' 



(0) 



A 



b 



3 

b 



2 3 6 a 



26 
3 



18) 



The leading order mass matrix is then diagonalized by 



in 



(0) 



^^TBMdiag (a + b, a, -a + b) U TBM . 



(19) 



Note that ruy is consistent with the TBM mass matrix form m™ M defined in Eq. (TTTT) since 
the identification A = a + |6, 5 = — 16 and D = b — a can be made. 

In order to modify the TBM mixing pattern, we introduce another flavon field 0, which 
transforms as 3 and couples to the lepton doublets via -^2X c <j) (£h u £h u ). Similar to £ and tpg, 
the unwanted couplings between right-handed charged leptons and <p are forbidden by the 
additional Z 3 symmetry. Different from the flavons ifx and <ps, the flavon field <fi is assumed 
to develop a VEV along the directions ((f)) = (0, — v^, v$). This vacuum alignment follows 
the orthogonality conditions ((f)) ■ (ipr) and ((f)) ■ ((/? ,<?) , where the "■" denotes the usual scalar 
product of 3- vectors. It has been shown in Ref. 22| that such an orthogonality condition can 
be realized within supersymmetry with "Lagrange multiplier" superfields, which are singlets 
under the flavor symmetry but couple to the flavon fields in the superpotential. The F-term 
conditions, which are equivalent to the orthogonality conditions, could then yield the desired 
vacuum alignments. 

The residual symmetry in the neutrino sector is now broken by (tps) — i v s, v s, v s) down 
to a Z2 symmetry G\. The additional term for m v from the extra flavon, after acquiring its 
VEV alignment ((p) = (0, —v^ v^), gives 



/o 



m 



(i) 



A 



-- 

2 w 



A 
2 * 



7 



(20) 



where c = 2x c ^. Writing C = \c, the TMi matrix structure given in Eq. (TTTT) is then 
reproduced. 

The above analysis could also be applied to other groups containing A4 as a subgroup. In 
the S4 model explored in Ref. 25j, the lepton doublets are assigned to a three dimensional 
representation 3i under S4 as well as two flavons ip ~ 2 and A ~ 3i. The vacuum alignments 
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are taken to be (ip) = (v^, ty,) and (A) = (v&, v&, v&). The neutrino mass matrix then reads 

(21) 



m 



(o) 



A 



/ 2 f d-f d-f\ 

d-fd + 2f -f 
\d-f -f d + 2f) 



where d = 2xdV^,/A and / = 2xfV&/A stem from the Yukawa coupling terms. This is 
equivalent to TBM, which can be seen by writing A = 2f, B = d — f and D — d + 3f. 
is diagonalized as 



Hence, 



m 



(o) 



^^TBMdiag i 



-d + 3/, 2d, d + 3/) Ur 



T 

TBM 



(22) 



Similar to the A 4 model, here we introduce a new flavon field (, which transforms as 3i 
but possesses a special vacuum structure, (£) ~ (0, — v^,v^), again possible to achieve with 
an orthogonality condition. The new flavon field leads to the following contribution to the 
neutrino mass term 



/o 



m 



(i) 



K 
A 



s 

s 2s 
-s 



-A 
o 

-2s J 



(23) 



where s = yv^/A with y being the Yukawa coupling between ( and the lepton doublets. 



Again, when rriv is added to , the TM4 mass matrix is reproduced. 

In these two examples, the VEVs of the new flavon fields should in principle be smaller 
than the other flavon VEVs, the reason being that C is proportional to sin 9 and expected 
to be suppressed with respect to the other entries in Eq. (flTT) . This could be achieved if the 
scale of new flavon fields is roughly one order of magnitude smaller than the others. 



,(°) 



IV. SUMMARY 



In this work we have studied an attractive neutrino mixing scheme TMi, in which the first 
column of the PMNS matrix has the same form as for tri-bimaximal mixing. The PMNS 
matrix can be described by using only two parameters: one rotation angle 9 and one CP 
phase ip. The physical observables, i.e. the three mixing angles and the Dirac phase, are 
therefore correlated via the two parameters, leaving us with rather definite phenomenology. 
While this was studied before, we noted here that the features that apparently emerge from 
global fits can be excellently described by this mixing scheme. Namely, if \U e s\ is non-zero, 
solar neutrino mixing is governed by sin 2 #12 = | — (|?7 e 3| 2 ), i.e. slightly less than |. If in 
addition the CP phase 5 is such that cos S is negative, and located around 7r, then atmospheric 
neutrino mixing is governed by sin 2 #23 = | — (|£/ e 3|), i.e. significantly less than |. These 
features, small negative deviations from |, large negative deviations from ^ and 5 around 
it, are (within la) the outcome of the global fit that we referred to. In fact, the sizable 
deviation from maximal mixing requires that the CP phase 5 is around tt, which can be 
tested in upcoming long-baseline neutrino facilities. We have also discussed potential flavor 
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symmetries behind the TMi scheme, and showed in particular that rather straightforward 
additions to existing models leading to tri-bimaximal mixing, be it A4 or S4, can lead to the 
mixing scheme. 

It will be interesting to see whether these features to which global fits seem to point survive 
the test of time. The mixing scheme that we studied here seems to be very well suited to 
describe the current data, and its rather simple structure adds to its attractiveness. 
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